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Abstract 



Elastic p photoproduction has been measured using the ZEUS detector at HERA. 
Untagged photoproduction events from ep interactions were used to measure the reaction 
IP ~^ P°P (p° ~^ 7r + 7r~) at photon-proton centre-of-mass energies between 60 and 80 GeV 
and \t\ < 0.5 GeV 2 , where t is the square of the four-momentum transferred at the 
proton vertex. The differential cross section da/dM nn , where M n7T is the invariant mass 
of the two pions, and the integrated cross section, a lp _, p o p , are presented; the latter was 
measured to be 14.7 db 0.4 (stat.) db 2.4 (syst.) ph. The differential cross section da/dt 
has an approximately exponential shape; a fit of the type A' t ex.p(— b' t \t\ + c' t t 2 ) yields a 
t-slope b' t = 9.9 dz 1.2 (stat.) dz 1.4 (syst.) fib. The results, when compared to low 
energy data, show a weak energy dependence of both a lp _, p o p and of the i-slope. The 
p° is produced predominantly with transverse polarisation, demonstrating that s-channel 
helicity conservation holds at these energies. 
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1 Introduction 



"Elastic" or "exclusive" production of p° mesons by photons, jp — > p°p, has been extensively 
studied in fixed target experiments up to photon-proton centre-of-mass energies W ~ 20 GeV, 
using both real and virtual photons |P|-||31||. Recently, the cross section for this reaction has also 
been obtained in an indirect measurement at the HERA ep collider, using quasi-real photons 
with space-like virtuality Q 2 between 4 • 10 -8 and 2 • 10~ 2 GeV 2 , at an average centre-of-mass 
energy (W) of 180 GeV |2|. 



At W values up to about 20 GeV, elastic photoproduction of p° mesons has the characteristic 
features of soft diffractive processes: the dependence of the production cross section on W 
is weak, the dependence on t, the square of the four-momentum transferred at the proton 
vertex, is approximately exponential, and the vector meson is observed to retain the helicity 
of the photon (s-channel helicity conservation, SCHC). Such energy and t dependences are 
also characteristic of hadronic diffractive processes. The similarity between p° photoproduction 
and hadronic processes can be understood in the framework of the Vector Dominance Model 
(VDM) ||33|| , in which the photon is assumed to fluctuate into a vector meson before interacting 
with the target nucleon; the reaction 7p — > p°p is thus related to the elastic process p°p — > p°p. 

At sufficiently high energies, diffractive interactions are usually described in terms of the ex- 
change of a pomeron, an object with the quantum numbers of the vacuum. Regge theory 
provides a framework in which many of the features of hadronic reactions can be described |34| . 
In particular, the energy dependence of diffractive processes is related to the intercept of the 
pomeron trajectory. Several models offer a description of diffractive vector meson produc- 



tion []35|]-[|I|; some of them are in the framework of perturbative QCD. The study of vector 
meson photoproduction at the energies available at HERA may thus help to clarify the nature 
of the pomeron. 

This paper reports a measurement of the elastic p° photoproduction cross section at (W) of 70 
GeV, based on about 6,000 ep — > ep7r + 7r~ events collected by the ZEUS experiment in 1993. In 
this measurement the final state electron and proton are not detected and the relevant kinematic 
quantities are determined from the measured three-momenta of the p° decay products, assuming 
that they are pions. 

The paper is organised as follows. After defining the variables relevant to p° production, we 
describe the experimental conditions and the event selection criteria, and then discuss the 
acceptance corrections and the background. From the analysis of the differential cross section 
da/dM wn , where is the invariant mass of the 7r + 7r~ system, we obtain the integrated 
cross section <y w ^ p o p . We then discuss the differential cross section da/dt and the angular 
distributions of the decay pions. Finally, from the value of da/dt at t = 0, the total p°p cross 
section is derived using the optical theorem and assuming VDM. 



2 Elastic p° photoproduction at HERA 

Elastic p° photoproduction was investigated by means of the reaction (see Fig. [I]) 

e(k) +p(P)^e(k')+p°(V)+p(P'), 



1 



where the symbols in parenthesis indicate the four-momenta of the particles involved. 

For unpolarised electrons and protons, two independent variables describe inclusive ep scatter- 
ing, since the ep centre-of-mass energy y/s = 2^J E e E p = 296 GeV is fixed by the energies of 
the electron (E e ) and of the proton (E p ) beams. The variables can be any two of the following 
four: 

• —Q 2 = q 2 = (k — k') 2 , the four-momentum squared carried by the virtual photon; 

• x = Q 2 /(2P ■ q), the Bjorken variable; 

• U = (<? ' P)/{k ' P) > the fraction of the electron energy transferred by the photon to the 
hadronic final state, measured in the proton rest frame; 

• W, the centre-of-mass energy of the 7*p system, where 

W 2 = (q + P) 2 = -Q 2 + 2y(k ■ P) + M 2 , 
M p being the proton mass. 

The hadronic final state, containing the scattered proton and the pions from the decay p° — ► 
7i + tc~ , is described by additional variables, including the invariant mass M„ of the two decay 
pions, the square of the four-momentum transferred at the proton vertex, t, and the polar and 
azimuthal angles, defined in section \IA[ of the decay pions in the nir centre-of-mass frame. 

For the data presented here, only the three-momenta of the final state pions were measured. 
Events in which the scattered electron was detected in the ZEUS calorimeter were rejected, 
thereby restricting Q 2 to be below Qmax ~ 4 GeV 2 . The median Q 2 is approximately 10~ 4 GeV 2 . 
To explain how the relevant kinematic quantities are obtained from the four-momenta of the 
two pions, we first consider the case Q 2 = Q^m = M 2 ^— (M e is the electron mass) and 
then discuss the effect of larger Q 2 . For Q 2 = Q^in (~ 10~ 9 GeV 2 for the kinematic range 
covered by the present data), the virtual photon is emitted with zero transverse momentum 
and with longitudinal momentum pz-y — —E 1 in the direction opposite to that of the proton 
beamF]. Energy and momentum conservation relate the photon energy E 1 to the two-pion 
system energy E nw and longitudinal momentum Pz-k-k by 2E^ ~ (E n7T — pz-nn) ■ The photon- 
proton centre-of-mass energy can then be expressed as: 

W 2 ~ AE^E P ~ 2[E^-pz^)E v . 

The p° transverse momentum squared in the laboratory frame, p^, approximates to —t: 

t = ( q - V ) 2 = -Q 2 -2q-V + M 2 n 
~ -2E^{E n7I + p z ^) + M 2 ^ 
- -(EL-Pz^) + M 2 n 

= ~P 2 T, 

1 Throughout this paper we use the standard ZEUS right-handed coordinate system, in which X = Y = Z = 
is the nominal interaction point, the positive Z-axis points in the direction of flight of the protons (referred to 
as the forward direction) and the X-axis is horizontal, pointing towards the centre of HERA. 
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where, in addition to 2E^ ~ [E^ — pz-Kw), we have used the approximation Q 2 = 0. Non-zero 
values of Q 2 cause p\ to differ from —t by ^Q 2 ; the effect on the measured distributions is 
discussed in section |5]. For this measurement, the minimum kinematically allowed value of \t\ 
is negligible, |t min | ~ 10~ 8 GeV 2 at W = 70 GeV. 

Fig. || shows the scatter plot of the reconstructed versus the true values of W and t for the sample 
of Monte Carlo events used to evaluate the acceptance (cf. section [5]). The difference of the 
reconstructed and the true value of W has a mean value of zero and an r.m.s. spread of 1.4 GeV. 
The analogous difference for t has also a mean value of approximately zero and an r.m.s. spread 
of 0.06 GeV 2 . The events far from the diagonal have Q 2 ^> . Throughout the analysis, 
the variable W was calculated using the above approximation. For t, the Q 2 dependence of the 
relation between t and p\ was taken into account in the acceptance correction, as discussed in 
section [5[ 

In the one photon exchange approximation, the ep and the 7*p cross sections for elastic p° 
production are related by 

'! + {!- y) 2 2(1 -y) Q\ 



d 2 o, 



ep~*epp 



,0 



a 



dydQ 



2ttQ 2 



y 



y 




7 p-*p"p 



y 



a 



(W,Q 2 )+ 



-p p 



(W,Q 2 ) 



(1) 



where a is the fine structure constant and dj. p 



p P (W, Q 2 ) and al p ~* p P (W, Q 2 ) are the re- 
spective cross sections for transversely and longitudinally polarised virtual photons. Following 



VDM, these are related by 



al p ~ 



\W,Q 2 



a 



TP- 



Q' 2 



(2) 



with 



a 



TP-*P P 



(W,Q 2 



a 



■yp^p 




^ / 1 + TH 



(3) 



where a lp ^ p o p {W) is the cross section for elastic photoproduction (Q 2 = 0) of p° mesons, and 



M p is the p meson mass. Substituting the latter two expressions into equation (ffl) yields: 

d 2 a P 



' ep^epp 

dydQ 2 



,0 



(4) 



where the function 



a 



2nQ 2 



l + (l-y) 2 2{l-y)(Q 



2 

min 



Q 2 



0^ 

Ml 



1 + 



(5) 



is the effective photon flux. 

The cross section a lp ^ p o p ((W)) for elastic p° photoproduction is thus obtained as the ratio of 
the corresponding acceptance corrected electron-proton cross-section, integrated over the y and 
Q 2 ranges covered by the measurement, and the effective photon flux Q 2 ) integrated over 
the same y and Q 2 ranges. This procedure determines the cross section, assuming VDM, for 
elastic production of p° mesons at Q 2 = averaged over the W range of the measurement. 
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3 Experimental conditions 



3.1 HERA 

During 1993, HERA operated at a proton energy of 820 GeV and an electron energy of 26.7 GeV; 
84 colliding electron-proton bunches were stored, with an additional 6 unpaired proton and 10 
unpaired electron bunches. These additional pilot bunches were used for background studies. 
The time between bunch crossings was 96 ns. Typical bunch currents were 10 mA for both the 
electron and the proton beam, providing luminosities of approximately 6 • 10 29 cm -2 s _1 . 



3.2 The ZEUS detector 



The components of the ZEUS detector are described in detail in ref. [0. A short description 
of those most relevant to the present analysis follows. 

Charged particles are tracked by the vertex detector (VXD) and the central tracking detector 
(CTD) which operate in a magnetic field of 1.43 T provided by a thin superconducting solenoid. 
The VXD |43| is a cylindrical drift chamber that surrounds the beam pipe and consists of 120 
radial cells, each with 12 sense wires. The VXD resolution in the plane transverse to the beam 
direction, for the data presented here, is 50 /im in the central region of a cell and 150 fim near 
the cell edges. The CTD consists of 72 cylindrical drift chamber layers, organised in 9 
superlayers covering the polar angle region 15° < < 164°. In 1993, the spatial resolution in 
the plane perpendicular to the beam was ~ 260 /mi. For the data presented in this paper, the 
combined CTD and VXD information provides resolutions for the primary ep interaction vertex 
of 1.1 cm in Z and 0.2 cm in the XY plane. The momentum resolution, for tracks traversing 
all superlayers, is cr(pt)/pt ~ \/(0.005) 2 p 2 + (0.016) 2 , where p t is in GeV. 



The high resolution uranium-scintillator calorimeter (CAL) |^5| consists of a forward (FCAL), 
a barrel (BCAL) and a rear (RCAL) part, respectively covering the polar regions 2.2° to 
36.7°, 36.7° to 129.1°, and 129.1° to 176.6°. The calorimeter parts are subdivided transversely 
into towers and longitudinally into electromagnetic (EMC) and hadronic (HAC) sections. A 
section of a tower is called a cell; each cell is viewed by two photomultiplier tubes. Holes of 
20 x 20 cm 2 in the centre of FCAL and RCAL accommodate the HERA beam pipe. From 



test beam data, the energy resolution was found to be <Je/E = 0.18/ J E(GeV) for electrons 



and <Je/E = 0.35/ y E(GeV) for hadrons. The performance, energy and time calibration of 
the calorimeter are continuously monitored using pedestal triggers, charge and light injection 
as well as the uranium radioactivity. The additional features relevant to the present analysis 
are the sub-nanosecond time resolution and the low noise of approximately 15 MeV for the 
electromagnetic and 25 MeV for the hadronic cells. 

The veto wall is used to tag events in which a proton has scattered off a residual gas molecule in 
the beam pipe ( "proton-gas" events) upstream of the ZEUS detector. It consists of two layers 
of scintillators, with overall dimensions of 500 cm x 600 cm, on both sides of an 87 cm thick 
iron wall centred at Z = —7.3 m. 

The C5 beam monitor, a small lead-scintillator counter assembly located at Z — —3.2 m, 
records the arrival times of halo particles associated with the proton and electron bunches 
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within 10 cm of the beam axis. It is used to verify the relative timing of the beams and to 
reject events due to proton-gas interactions. 

The luminosity detector (LUMI) |Q measures the luminosity by means of the Bethe-Heitler re- 



action ep — > ejry; a detailed description of the method used is given in |32fl . The bremsstrahlung 
events are identified by measuring the radiated photon in a lead-scintillator calorimeter placed 
in the HERA tunnel downstream of the interaction point in the direction of the outgoing 
electron. 



3.3 Untagged p photoproduction trigger 



ZEUS uses a three level trigger system | 42fl . The data presented here come from the "untagged 
photoproduction trigger", designed to select vector meson photoproduction events [57]. The 
term "untagged" refers to the fact that the scattered electron escapes undetected through the 
beam pipe hole in the RCAL and is not detected in the LUMI detector. 



The trigger conditions can be summarised as follows: 



• First-level-trigger (FLT): 

— At least 464 MeV deposited in the electromagnetic section of RCAL. 

— At least one track candidate in the CTD. 

— Less than 3750 MeV deposited in the calorimeter towers surrounding the beam pipe 
in the forward direction. This requirement suppressed proton beam-gas events and 
a significant fraction of the photoproduction cross section. 

The trigger was vetoed if hits were present in the C5 or in the veto wall counters, with 
timing consistent with that of a p-gas collision occurred upstream of the interaction point. 

The resulting FLT rate was ~ 10 Hz at a luminosity of 6 • 10 29 cm -2 s -1 . 

• Second-level-trigger (SLT): 

— Events with calorimeter timing indicating that the interaction had occurred up- 
stream of the interaction point were rejected. 

• Third-level-trigger (TLT): 

— Cosmic ray events were discarded on the basis of calorimeter timing. 

— A tighter calorimeter timing rejection was applied. 

— A cut on the Z value of the reconstructed event vertex of ± 85 cm was imposed. 

The rate of events passing the untagged photoproduction trigger at the third level was 
about 1.2 Hz. For some fraction of the data-taking period a factor of 3 prescale was 
applied. 



The requirements that a signal be detected in RCAL and a track be seen in the CTD effectively 
selected events with photon energies between 0.5 and 4 GeV, corresponding to 40 < W < 120 
GeV. 
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4 Event selection 



During 1993, approximately 7 ■ 10 6 events were recorded, corresponding to a total integrated 
luminosity of about 550 nb _1 . The data presented in this paper correspond to an effective 
luminosity, accounting for the effects of the prescaling of the trigger mentioned in the previous 
section, of (240 ± 8) nb" 1 . 

The following offline requirements were imposed in order to obtain the final p° photoproduction 
sample: 

• exactly two tracks from particles of opposite charge and both associated with a recon- 
structed vertex; 

• transverse momentum greater than 200 MeV and hits in at least the 3 innermost CTD 
superlayers for each of the two tracks, thus restricting the data to a region of well under- 
stood track reconstruction efficiency; this restriction approximately translates into one on 
the track pseudorapidity (rj = — In tan 0/2) of \r]\ ^ 1.8; 

• vertex position within ± 40 cm of the nominal interaction point and within a radial 
distance of 1.5 cm from the beam axis (the interaction region was centred at Z = —6 cm 
and had an r.m.s. width of ±11 cm); 

• total energy in the forward calorimeter Efqal < 1 GeV, thereby limiting the contamina- 
tion from the reaction ep — > ep°X, where X is a hadronic state of mass Mx into which 
the proton had dissociated; 

• no energy deposits larger than 200 MeV in any calorimeter cell outside a circular region 
around the track impact point with a radius of 30 cm in the EMC and 50 cm in the HAC, 
thus rejecting events with neutral particles or with charged particles outside the region 
of sensitivity of the tracking system. 

A total of 13570 events satisfied the above criteria. 

The pion mass was assigned to each track and the analysis was restricted to events reconstructed 
in the kinematic region defined by: 

60 < W < 80 GeV, 
0.55 < M n7T < 1.0 GeV, (6) 
p\ < 0.5 GeV 2 . 

In the chosen energy range the acceptance is well understood. The restricted mass range reduces 
the contamination from reactions involving other vector mesons, in particular from elastic </> 
and to production, as well as from photon conversions. The restricted px range reduces the 
contamination from events with diffractive dissociation of the proton (ep —>■ ep°X). The final 
sample contains 6381 events. 

The invariant mass spectrum is shown in Fig. |3] before and after the offline selection. The data 
are dominated by the p° signal. The corresponding mass spectra of like-sign two track events 
are also shown as the shaded areas. The small peak just above the tttt threshold in Figs. |3|b 
and c is due to — > K + K~ events, where the pion mass has been erroneously assigned to the 
tracks. 
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5 Monte Carlo generators and acceptance calculation 



The reaction ep — > ep°p was modelled using two different Monte Carlo generators. 



The first generator, DIPSI |48|], describes p° photoproduction in terms of pomeron exchange. 



Based on the model of Ryskin P?|, it assumes that the exchanged photon fluctuates into a qq 



pair which then interacts with the pomeron emitted by the incident proton. The pomeron is 
described in terms of a gluon ladder. The cross section is proportional to [a s (g 2 )] 2 ■ [xg(x, q 2 )] 2 , 
where a s (q 2 ) is the strong coupling constant, xg(x,q 2 ) is the gluon momentum density in 
the proton, x is the fraction of the proton's momentum carried by the gluon ladder and 2q 2 , 
in the leading logarithm approximation, is the upper limit for the virtuality of the two t- 
channel gluons of the gluon ladder. Once a s and the gluon momentum density are fixed, 
the W and t dependences are determined. The process under study is sensitive to values of 
x « Mp/W 2 m 1CT 4 and q 2 m 0.25 M 2 « 0.15 GeV 2 . The latter is below the expected region of 
validity of the calculation; a parametrisation for the product [a s (g 2 )] 2 • [xg(x, q 2 )] 2 was however 
found for which the model describes all measured distributions well [H9J . For this choice of the 
parameters the cross section has a very weak dependence on W. The two-pion invariant mass 
M n7T was generated so as to reproduce, after reconstruction, the measured distribution. 



The second generator, LEVEC, was developed within the HERWIG framework ||50|| . It assumes 
expression (|3|) for the Q 2 dependence of the cross section; the contribution of longitudinal 
photons is neglected. The generated events were weighted such that all other distributions 
(i.e. those over W, M n7T , p\ and the angular distributions of the decay pions), after detector 
simulation and event reconstruction, have the same shape as those of the data. 

For both programs, the angular distribution of the decay pions was assumed to be that expected 
on the basis of SCHC Elf . 



The acceptance for elastic p° production was calculated using both the DIPSI and the LEVEC 
generators. Fig. [|a, b and c respectively show the acceptance as a function of M nn , W and p\. 
The acceptance includes the effects of the geometric acceptance of the apparatus, of the detector 
efficiency and resolution, and of the trigger and reconstruction efficiencies. The trigger efficiency 
is ~ 43%. The average acceptance is about 7%. The acceptance increases with increasing tttt 
mass, has a broad maximum for W ^ 70 GeV and is almost independent of the p° transverse 
momentum squared. In order to convert the measured dN/dp 2 ^ distribution to the differential 
cross section da/dt, the p\ acceptance (Fig. |]c) was multiplied by a correction factor F, which 
is, bin by bin, the ratio of the p\ and t distributions at the generator level. Fig. [|d shows F as 
a function of p\. 

To produce the acceptance corrected cross sections, the generator DIPSI was used. The differ- 
ence between the results for the acceptance obtained with the two generators was taken as an 
estimate of the systematic error due to the model dependence of the acceptance calculation. 

From the comparison of the distributions of the reconstructed and generated events, the 
resolutions in and p\ were found to be cr Mjrjr ~ 20 MeV - consistent with the mass 



widths found for neutral strange particles reconstructed using the CTD information [J52] - and 



a p 2 T ~ 0.01 GeV 2 . 

The DIPSI and LEVEC simulations were also used to produce samples of elastic u and 
photoproduction events for the study of the contamination from such processes. 
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The reaction ep — > ep°X, where X is the hadronic state resulting from the diffractive dissoci- 
ation of the proton (see section was simulated using PYTHIA | 53fl . A cross section of the 
form d 2 'a / ' dtdM\ oc e bt /M x , with b = 4.5 GeV -2 , was assumed; the maximum value of Mx 



was fixed by the condition M x /W 2 < 0.1 ||54|| . The exponent (3 was varied between 2 and 3, 
consistent with the result f3 = 2.20 ± 0.03 recently obtained at Fermilab |55[ for the diffractive 
dissociation of the proton in pp collisions. A second generator (RHODI) was also used, based 
on the model calculation of ref. |56||. In both generators the p° decay angular distributions were 



assumed to be the same as those of the elastic events. 



6 Backgrounds 

After applying the selection criteria described in section |], the data still contain contributions 
from various background processes to the reaction ep — > e7r + 7r~p: 

• inelastic p° production, ep —>■ ep°X, in which the proton diffractively dissociates into a 
hadronic state X not detected in the calorimeter; 

• elastic production of u> and (f) mesons; 

• beam-gas interactions. 

In order to estimate the contamination from the inelastic channel ep ep°X, the energy 
distribution in the forward calorimeter (E FC al) was studied. Fig. [5] shows the distributions 
for both elastic and inelastic events as obtained from Monte Carlo simulations based on the 
generators described above, along with the distribution for the data. To obtain these plots, 
the cut Epcal < 1 GeV was not applied. The plot for elastic Monte Carlo events (Fig. ||a) 
shows the FCAL energy spectrum resulting from noise (mainly the uranium radioactivity). 
Nearly all events have Epcal < 1 GeV. For the inelastic Monte Carlo sample as well as for 
the data, the Epcal spectrum extends to much higher values. Energy deposits in the forward 
calorimeter greater than 1 GeV were therefore ascribed to inelastic reactions in which part of 
the diffractively produced hadronic state X was detected in FCAL. The number of residual 
inelastic events in the data with FCAL energy smaller than 1 GeV was then estimated as 

The number of Monte Carlo and data events, {N [Epcal > 1 GeV)}A/c an d {N (Epcal > 
1 GeV)} data, were computed in the region 1 < Epcal < 8 GeV, where the inelastic Monte 
Carlo describes the data well. The overall contamination integrated up to \t\ = 0.5 GeV 2 was 
estimated to be 11% ±1% (stat.) ±6% (syst.). This was obtained using PYTHIA with (3 = 2.5. 
The systematic error reflects the sensitivity of the result to the value of the exponent f3, which 
was varied between 2 and 3, and to the use of RHODI instead of PYTHIA. As a check, the 
sum of the elastic and inelastic Monte Carlo distributions of Fig. |5|a and [|b were fitted to 
the data of Fig. ||c; the normalisations of the simulated distributions were free parameters of 
the fit. The x 2 °f the fit has a broad minimum around (3 = 2.5. This method gave a 12% 
contamination, consistent with the result given above. The result depends very little on the 
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shape of the generated t distribution: the estimate of the contamination varies by less than the 
quoted statistical error for a change of the exponential slope b between 4 and 6 GeV -2 . 

The contamination was also studied as a function of t and W. It was found to vary from 6% 
for \t\ < 0.05 GeV 2 to 19% for \t\ ~ 0.5 GeV 2 ; it increases by 3 ± 2% as W increases from 65 
to 75 GeV. 

The contamination due to the elastic production of <fi and u mesons was estimated from Monte 
Carlo simulations. The selection cuts described in section |] were applied to the simulated events 
after reconstruction. As an example, the contamination due to us production was estimated as 

A^o^ 
A^c^ + ApUp 

where is the acceptance for elastic lu events. Assuming a cross section ratio of o w jo p = 
0.1 [p7| , a contamination of (1.3 ± 0.2)% is obtained. A similar procedure applied to events 
results in an estimated contamination of (0.3 ± 0.1)%, mainly due to the — ► 3n decay mode. 
These contributions were not subtracted but were included in the systematic error. The con- 
tributions from inelastic u and production were negligible. 

Electron beam-gas and proton beam-gas contaminations were deduced from the pilot bunch 
event samples to which the cuts described in section ^ were applied. The number of events 
passing the cuts was then scaled by the ratio between the electron (proton) current in the 
paired bunches and the current in the electron (proton) pilot bunches. The contamination due 
to electron-gas interactions was estimated to be 2.3 ± 0.5%, while that due to proton-gas events 
was found to be 0.3 ± 0.2%. 

All subsequent results are shown after subtraction of the contributions from inelastic proton 
diffraction and beam-gas interactions. 



7 Results 

7.1 Differential cross section da/dM n7r 

Fig. ^| shows the acceptance corrected differential cross section da/dM nn . The mass distribution 
is skewed compared to a Breit-Wigner distribution: there is an enhancement of the low mass 
side and a suppression of the high mass side. This distribution can be understood in terms of 
the interference between the resonant 7r + 7r~ production and a non-resonant Drell-type back- 
ground |BB[ as discussed by Soding ]59]|. In order to extract the contribution of the resonant 



part of the differential cross section da/dM nn , we followed a procedure similar to that described 



in refs. || [Tjj, p5|| . The function 

da/dM^ = f p ■ BW P {M^) + fj ■ I(M n7T ) + f PS (7) 
was fitted to the measured mass distribution. The term 

BW(M \ M 7T7T MpTp(M n7T ) 

Pl ™ j - _ + M 2 T 2 {Mn7T) W 
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is a relativistic Breit-Wigner function, with a momentum dependent width |5D 




W = r r- -JL, (9) 



where To is the width of the p°, p* is the it momentum in the irn rest frame and p* is the value 
of p* at the p° nominal mass M p . The function 

M 2 - M 2 

KM ) = p - — (10) 

1 ™ } {ML - M 2 ) 2 + M 2 p Y 2 p {M^) 1 ) 

is a parametrisation of the interference term. The background term fpg was taken to be 
constant. The free parameters in the fit were M p , T and the coefficients f p , fi and fpg. 

The results of the fit are presented in table H and in Fig. |^. The x 2 /ndf is 1.4, for ndf— 13. 
The fitted values of the p° mass and width are in good agreement with the accepted ones |6T] . 
The background term fpg is consistent with zero, within a large error; similar results for fpg 
were obtained by earlier experiments |], |1(|| using the functional form (]?]). 



The contribution of the resonant term increases with \t\, ranging from 86% of the events for 
|*| = 0.01 GeV 2 to 95% for \t\ = 0.5 GeV 2 . 

The interference and the background terms were also studied as a function of W and of the 



decay pions' angular variables, cos 6^ and <fih, defined in section |7l| No dependence on these 
variables was found. 

The fit was repeated using different assumptions for the functional form of da/dM nn . 

• Parametrisation (0) is only an effective one, as it leaves the interference term independent 
of the resonant and non-resonant terms, which is strictly speaking inconsistent with the 
Soding mechanism. We therefore fitted the following functional form to the invariant 
mass distribution: 



da/dM^ 



A 



M^M P Y P 



M 2 W - M 2 + iMpYp 



+ B 



(11) 



where A, B, M p and Tq were free parameters of the fit. For V p expression (||) was used. 
The non-resonant amplitude B was taken to be constant and real; it was also constrained 
to be non-negative. The results for the parameters are given in table 0. The x 2 /ndf was 
1.4, with ndf = 14. 

The following alternative expressions for the width of the p° were adopted in the func- 
tions fBIJIT 



T p (M n 



(12) 



r P (M w 



r„ 



v _ 



i + (p*/p*o¥ 



(13) 
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The Breit-Wigner was parametrised, following refs. p"8l B(J, as 



nw(M \ - 1 M nw M p T p (M nn ) 

p[ ™ } ~ P* (M^ - Mjy + M^ p (M n7T ) li4j 

and expression (|13|) was used for the width. 
The parametrisation given in ref. was used: 

dxr/dM m = f p ■ BW p (M nn ) • (l + d \(M p /M n7T ) 2 - l] + C 2 \(M p /M n7T ) 2 -if). (15) 



The fit was repeated for the three mass dependent widths (|9]), (|12"D and (|T3|). 
• The phenomenological parametrisation proposed by Ross and Stodolsky |62| was used: 
da/dM nn = f p ■ BW p (M n7r ) ■ (M p /M nn ) n + f PS , (16) 

where the factor (M p /M nn ) n accounts for the skewing of the shape of the p° signal. The 
term fp$ was taken to be constant. Here again the fit was repeated for the three mass 
dependent widths (|S]), ( |T2"D and ([□]). The parameter n was found to be n = 4.9 ± 0.5, 
n = 5.8 ± 0.5 and n = 4.9 ± 0.5 for the three forms of the width, respectively. 

In none of these cases did the quality of the fit change appreciably, as can be seen from table |^, 
in which the values of x 2 /ndf obtained for the various functional forms are summarised. The 
fitted values of the p° mass and width varied from 763 to 772 MeV and from 141 to 155 MeV, 
respectively. As we discuss in the next section, the values of the resonant part of the cross 
section were quite stable. 



7.2 Integrated — > p°p cross section 

The cross section a lp ^, w + n - p at Q 2 = 0, integrated over the M nn and t regions specified in 
and averaged over the range 60 < W < 80 GeV, can be obtained from the data as 

-A-tt + tt- — 



'7£> — >7T~TT p 



where N n + n - is the number of observed events after background subtraction, e is the overall 
acceptance, L is the effective luminosity and $ = 0.02419 is the effective photon flux factor 
(see eq. |]) integrated over the specified W and Q 2 ranges. In order to extract the cross section 
for the resonant process 7^ — > p°p, it was assumed that the p° meson decays to 7r + 7r~ with a 
100% branching ratio and the fit procedure described in section |7.1| (with expressions (|7|-[ToD) 
was used. The resonant part of the total 7r + 7r~ signal is given by the parameter f p multiplied 
by the integral of the relativistic Breit-Wigner curve, that is the area under the dotted curve 
in Fig. There is some arbitrariness in the choice of the integration limits of the Breit- 
Wigner curve. The integral was carried out in the range 2M n < < M p + 5r , where the 
p° mass and width values were taken from the fit; the quantity M n is the pion mass. This 
requires an extrapolation beyond the measured region. The upper limit for the integration 
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range approximately corresponds to the mass of the nearest resonance, the p(1450), with the 
same quantum numbers and quark content as the p°. The value of the cross section, for 
2Af w < Mm < M p + 5r , \t\ < 0.5 GeV 2 and averaged over the region 60 < W < 80 GeV, was 
measured to be: 



a 7 p_ p o p = 14.7 ± 0.4 (stat.) ± 2.4 (syst.) fib. (17) 

The systematic error is dominated by the uncertainties on the acceptance determination (13%), 
on the number of p° signal events (7%), and on the inelastic background determination (7%). 
If the integration is limited to the measured region, 0.55 < M^ < 1 GeV, the cross section 
value is 12.4 fib, i.e. lower by a factor £ = 1.19. If the integral is computed up to M p + 4T , the 
cross section decreases by 3%; if instead it is extended to M p + 6r , the cross section increases 
by 2%. 

The uncertainty on the acceptance determination has three main contributions: 



1. the uncertainty on the calorimeter trigger efficiency near the threshold (9%); 

2. the difference between the results obtained with DIPSI and with LEVEC (8%); 

3. the sensitivity of the results to the cuts, notably that on the minimum number of CTD 
superlayers traversed by each track (6%). 



The various alternative functional forms for da/dM nn described in the previous section were 
also used to extract the resonant part of the signal. The values obtained were centred around 
the result given in ([TJ]) but spanned the range a 1P ^ p o p = 13.6-15.4 fib, corresponding to a 

3%. The method of Spital and 



1|% maximum variation. The corresponding variation of £ is 
was also used to obtain the cross section, yielding 



Yennie 63 



a jp-*p°p 2 



7rr da 



dM^ 



15.5 ±0.4 fib, 



Mo 



(18) 



where the p° mass and width were those given in table [TJ. The result obtained with the Spital 
and Yennie method depends linearly on the value used for To; it is also sensitive to the value 
of M p , since da/dM n7T is a steep function of M n7T in the region M 7r7r ^750 MeV, as can be seen 



from Fig. y. If the values of M p and Tq given in ref. |61j are used, the cross section changes by 
less than 1%. If T is kept fixed and M p is varied between 760 and 780 MeV, the corresponding 
change in the cross section is 23%. 

Table |3| summarises the results. We have taken the spread into account by including a ± 7% 
contribution to the systematic uncertainty of the cross section. 

The effect of real photon radiation by the incoming or the outgoing electron and that of vacuum 
polarisation loops is to lower the measured value of the cross section. The size of the correction 



was estimated to be smaller than 4% [64]. The correction was not applied; instead a 



contribution was added to the systematic uncertainty. 

Table § summarises the contributions to the systematic error. The total systematic error was 
obtained by summing all contributions in quadrature. 



12 



Result (|TTD for the cross section a lp ^ p o p is presented in Fig. 0a together with a partial compi- 
lation of low energy measurements. The figure only includes results explicitly corrected for the 
interference term and the non-resonant background. We also show the ZEUS result f3"2" |, ob- 
tained indirectly using tagged photoproduction data from the 1992 data-taking period. Fig. |7|b 
shows the result ( |T8| ) obtained with the method of Spital and Yennie, along with a compilation 
of low energy results obtained with the same technique. The dashed curve in both figures is 
a parametrisation by Schuler and Sjostrand |55| ; the W dependence of the curve is based on 
Regge theory and on the assumption that the intercept of the pomeron is a(0) = 1 + 0.0808 
(the soft pomeron). The same general trend of the data is seen in the two figures. There 
are however differences in the results obtained by individual experiments; these differences at 
least in part reflect the ambiguity in the definition of the p° production cross section due to 
the finite width of the p°. The comparison between different experiments - and between the 
experimental results and the theoretical expectations - should thus be taken with caution. The 
curve satisfactorily reproduces the energy dependence of the data. 



7.3 Differential cross section da/dt 

Fig. |8]a shows the acceptance corrected differential cross section da/dt, integrated over the 
p° mass region 2M n < M n7T < M p + 5Fo. It was obtained by multiplying the differential 
cross section for the measured range 0.55 < M n7T < 1 GeV by the factor £ = 1.19 discussed 



in section [T.2| . This assumes that, in each t bin, the ratio of the integral of the relativistic 
Breit-Wigner distribution over the range 2M n < M^ n < M p + 5r to that over the range 
0.55 < M,,-,!- < 1 GeV is the same, i.e. that the mass and the width of the p° are the same in each 
bin. The contamination from inelastic p° production was subtracted bin by bin. Background 
and interference terms were also subtracted; their contribution was found by repeating the 
fitting procedure described in section |7.1| in each t bin, fixing the values of the p° mass and 



width to those given in table [I]. Fig- H shows the results obtained for da/dt by applying the 
Spital and Yennie method in each t bin. 

The data were fitted with the function 



wr A '- e ' bM (19) 



in the range \t\ < 0.15 GeV 2 and with the function 



da 
d\t\ 



A' t -e 



(20) 



in the range \t\ < 0.5 GeV 2 . Both functions describe the data well. The results of the fits using 
expression Q2"0| ) are shown in Fig. |8|. Table |5] gives the values of the parameters obtained in the 
fits. The difference between the results of the fits to the points of Fig. |]a and b was taken as 
an estimate of the systematic uncertainty on the determination of the resonant fraction of the 
cross section in each bin. The other contributions to the systematic errors on A t , A' t and bt, 
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b' t are the uncertainties on the acceptance and on the inelastic background determination. All 
contributions to the systematic errors were summed in quadrature. If the first bin is excluded 
from the fit to the points of Fig. |]a, the values of the slopes b t and b' t increase by 9% and 7%, 
respectively; the variation is smaller for the fit of Fig. |8|b. 

Fig. H shows the result for the slope b' t as obtained from fitting eq. ( p0|) to the points of 
Fig. |^a together with a compilation of fixed target results. All the results shown were obtained 
from fits of the form (f2"0"|); the data were explicitly corrected for interference and non- resonant 



background, with the exception of those of ref. ||. The data of refs. || |T!| have the somewhat 
large minimum \t\ values of 0.2 and 0.15 GeV 2 , respectively. For W > 4 GeV, the results 
do not depend strongly on W, as expected from Regge theory, which predicts a logarithmic 
dependence of the slope on W if one trajectory dominates. 



7.4 Decay angular distributions 

The p° decay angular distributions can be used to determine elements of the p° spin-density 
matrix 115111. 



In the s-channel helicity frame, in which the p° direction in the photon-proton centre-of-mass 
frame is taken as the quantisation axis, the decay angular distribution H (cos 9h,4>h,&h) is a 
function of the polar angle 9h of the 7r + in the p° centre-of-mass frame, of the azimuthal angle 
4>h between the decay plane and the 7-p plane (the p° production plane) and of the angle 
between the p° production plane and the electron scattering one. For t = t min , the photon and 
the p° are collinear and <ph is not defined. 

In the present experiment, the lepton scattering plane is not measured since neither the recoil 
proton nor the scattered electron are detected. Furthermore, the azimuthal angle <ph can be 
determined only if the direction of the virtual photon is approximated by that of the incoming 
electron. It has been verified by Monte Carlo calculations that this is a good approximation. 
The experimental resolution in 0, is approximately 40°; it is a function of t and improves with 
increasing \t\. The resolution in cos 9h is ~ 0.03. In the following we present the results for 
the one dimensional distributions H cos6h (cos9h) and obtained from H(cos9h, (fih, &h) 

after integrating over 0,, $^ and over cos 6^, respectively. For unpolarised or transversely 
polarised electrons and a J p = 1~ state decaying into two pions, the functions H cos9h (cos9 h ) 



and H^ h (cf)h) can be written as [^T], p| : 



H co *e h = \[l~ r ° 4 + (3r ° 4 - 1) cos 2 9 h ], (21) 

^ = ^(1-2^ cos 20,), (22) 

where Tqq and r°i x are p° density matrix elements. In particular gives the probability that 
the p° is longitudinally polarised. Assuming s-channel helicity conservation (SCHC), Tqq can 
be related to R, the ratio of the p° production cross sections for longitudinally and transversely 
polarised virtual photons |2Tf : 

1 r 04 

^=-7^04' (23) 
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where e is the virtual photon polarisation, i.e. the ratio of the longitudinally to transversely 
polarised photon fluxes. The present data have e = 0.998, essentially constant over the W 
range covered by the measurement. 



Fig. [K] shows the differential cross sections da/dcosOh and da /d4>h- The curves are the results 
of the fits of the functions (21) and (22) to the data. A dominant sin 2 Oh contribution is visible, 



indicating that the p mesons are mostly transversely polarised. This is reflected by the result 

,,04 
00 



of the fit, which yields Tqq = 0.055 ± 0.028, where the error is statistical only. The fitted value 



of 

r i—i is 0.008 i 0.014, consistent with zero, as expected on the basis of SCHC. If one uses 
expression ( p3|) to determine R from the fitted value of r^, one obtains R = 0.06 ± 0.03. The 
average value of Q 2 for the data discussed in this paper, computed assuming the Q 2 dependence 
given in (|j), is ~ 0.1M 2 = 0.05 GeV 2 . This gives, using expression (fj), R = 0.1, consistent 
with our result. 



7.5 Total p°p cross section 



By using the vector dominance model and the optical theorem, the measured value of da /dt at 
t = can be used to obtain the p°p total cross section. 



For p°p scattering the optical theorem reads: 



dt 



l + rf 
16tt 



(24) 



where r] is the ratio of the real to the imaginary part of the forward p°p scattering amplitude. 
We assume that rj = 0. On the other hand, within VDM, elastic p° photoproduction is related 
to the elastic p°p cross section; in particular, for t = 0: 



da^p^pQp 



dt 



Ana da p 0p^ p o p 



f 2 



dt 



(25) 



t=o 



where Ana / f 2 is the probability for the 7 
p. 393). 



p° transition. We take f 2 /An = 2.20 (cf. e.g. B 



Using the intercept da lp . 
equations (JT~ 



*p° P /dt\ 



t=o 



133±11±27 pb/GeV 2 given in section |7!3] , and combining 
and (IM), one finds a to t{p°p) = 28.0 ± 1.2 (stat.) ± 2.8 (syst.) mb, where the 



errors reflect only the uncertainty on the measured value of the intercept. The systematic error 
does not include the uncertainties from the model dependence of the assumptions made nor 
from the values of 77 and Air/ f 2 . 

The result is consistent with those found at lower energies (see e.g. [[j], 0])- It is also in 
agreement with the expected value of 27.8 mb at W = 70 GeV, obtained from the parametri- 

: 13.63(VU 2 ) - 0808 + 31.79(iy 2 )- - 4525 , based on 
and on fits to Tip data at centre-of-mass 



sation |5| of a$ p w l/2[a tot {n + p) + a tot (ir~p)} 
the soft pomeron, the additive quark model |6 
energies ranging between 6 and 25 GeV. 
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8 Summary 



The ZEUS detector at HERA has been used to study the photoproduction process 7p — > p°p. 
The integrated cross section as well as the differential cross sections da/dM n7T and da/dt at an 
average photon-proton centre-of-mass energy (W) = 70 GeV have been measured. 

The 7T7T mass spectrum is skewed, compared to a relativistic Breit-Wigner distribution, as also 
observed at low energy. The integrated 7p — > p°p cross section, for 60 < W < 80 GeV, 
|t| < 0.5 GeV 2 and 2M n < M n7T < M p + 5r , is 14.7 ± 0.4 (stat.) ± 2.4 (syst.) ph. This 
result, in conjunction with the measurements at low energy, is consistent with the weak energy 
dependence expected on the basis of Regge theory and a pomeron intercept of 1.08 (the soft 
pomeron). This is at variance with the behaviour of the cross section for elastic photoproduction 
of J /ip mesons [6J| and elastic production of p° mesons for Q 2 > 7 GeV 2 J70J. 



The differential cross section da/dt has an approximately exponential shape with a slope con- 
sistent with the results obtained by low energy experiments with W > 4 GeV. This is also in 
accord with Regge theory, which expects a logarithmic dependence of the slope on W. 

The p° decay angular distributions have been studied. The p° mesons are mainly produced in a 
transversely polarised state; the probability to find the p° with longitudinal polarisation is Tqq = 
0.055 ±0.028. The value of fjt.i is 0.008 ±0.014, in accord with s-channel helicity conservation 
(SCHC). If SCHC is assumed, one obtains R = 0.06 ± 0.03 for the ratio of the p° production 
cross sections for longitudinally and transversely polarised virtual photons, consistent with the 
value expected if VDM is assumed. 

From the measured value of da/dt at t = 0, within the VDM framework and by using the optical 
theorem, a value 0"*°^ = 28.0 ± 1.2 (stat.) ± 2.8 (syst.) mb for the total p°p cross section at 
W = 70 GeV is found, in agreement with extrapolations of the pion-proton total cross section 
data based on Regge theory and a soft pomeron. 
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Parameter value stat. error 



M p 


0.764 GeV 


0.003 GeV 


r 


0.142 GeV 


0.007 GeV 


fp 


9.86 fib 


0.23 fib 


h 


3.47 fib GeV 


0.34 fib GeV 


fps 


0.0 fib/GeV 


5.5 fib/GeV 



Table 1: Results of the fit to the mass spectrum for 60 < W < 80 GeV and |t| < 0.5 GeV 2 
(expressions ([7|-|l(])). Only statistical errors are given. 



Parameter value stat. error 

0.764 GeV 0.003 GeV 

0.148 GeV 0.007 GeV 

-3.10 fib 1 / 2 0.04 fib 1 / 2 

2.0 fib 1 / 2 GeV- 1 / 2 0.2 fib 1 / 2 GeV- 1 / 2 



Table 2: Results of the fit to the mass spectrum for 60 < W < 80 GeV and \t\ < 0.5 GeV 2 
(expression flTTD). Only statistical errors are given. 



M p 

r 

A 
B 
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Functional form 




XVndf 


(ndf) 


0.55 < M nw < 1 GeV 


cr^p^pOp [fxb] 
2M V < M n7r < M p + 5r 


Expressions 


(EH 


,00) 


1.46 


13) 


12.4 ± 0.3 


14.7 ± 0.4 


Expressions 






1.64 


;i3) 


11.5 ± 0.2 


14.5 ±0.3 


Expressions 


(Hi 


,0,0) 


1.45 


;i3) 


12.3 ± 0.6 


14.3 ±0.7 


Expressions 


(hi 


1,0,0) 


1.46 


;i3) 


12.1 ± 0.6 


13.6 ±0.6 


Expressions 


(OB) 


1.38 


;i4) 


11.9 ± 0.3 


14.3 ± 0.4 


Expressions 


flUi 


1) 


1.45 


13) 


12.4 ±0.5 


14.7 ±0.6 


Expressions 


Hi 


m 


i /in 

1.49 


'1 Q\ 

,13) 


ion _i_ n c 

12.0 ± 0.5 


15.4 ± 0.6 


Expressions 


m 


m 


1.47 


13) 


12.3 ±0.5 


14.4 ±0.6 


Expressions 


HI) 


1.45 


13) 


12.3 ±0.6 


14.6 ±0.8 


Expressions 


(£1 


0) 


1.45 


13) 


11.9 ±0.6 


15.2 ±0.8 


Expressions 


H, 


13) 


1.46 


13) 


12.3 ±0.6 


14.3 ±0.7 


Expression I 


PD 










15.5 ±0.4 



Table 3: Values of x 2 /ndf and results for the cross section 7p — > p°p obtained by fitting the 
different functional forms described in the text to the data of fig. [f| The indicated errors are 
the statistical ones only. 



Contribution from Error 

Luminosity 3% 

Acceptance: trigger efficiency 9% 

Acceptance: model dependence 8% 

Acceptance: sensitivity to cuts 6% 

Inelastic background subtraction 7% 

Background due to elastic u and <fi production 1% 
Procedure to extract the resonant part of the cross section 7% 

Radiative corrections 4% 

Total 17% 



Table 4: Individual contributions to the systematic error on the integrated cross section. 
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Data from Fig. |a 


Data from Fig. |b 


A t 


= 139 ± 6 (stat.) ± 26 (syst.) 


/xb/GeV 2 


A t 


= 159 ±6 (stat.) ±28 (syst.) 


Atb/GeV 2 


b t 


= 10.4 ± 0.6 (stat.) ± 1.1 (syst.) 


GeV- 2 


h 


= 11.2 ± 0.5 (stat.) ± 1.1 (syst.) 


GeV" 2 


A' t 


= 133 ± 11 (stat.) ± 27 (syst.) 


Aib/GeV 2 


A't 


= 155 ± 14 (stat.) ± 29 (syst.) 


/ib/GeV 2 


K 


= 9.9 ± 1.2 (stat.) ± 1.4 (syst.) 


GeV- 2 


b' t 


= 11.1 ± 1.5 (stat.) ± 1.4 (syst.) 


GeV" 2 


< 


= 2.3 ± 2.8 (stat.) 


GeV- 4 


c't 


= 4.7 ±3.3 (stat.) 


GeV- 4 



Table 5: The results of the fits to the data of Fig. |8|a (left) and Fig. |8|b (right) with function (|20|). 
For an explanation of the symbols see text. 
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Figure 1: Elastic p production in ep collision; the exchanged virtual photon is indicated as 7*. 
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55 60 65 70 75 80 85 0.1 0.2 0.3 0.4 0.5 0.6 



W gen (GeV) llt g J (GeV 2 ) 

Figure 2: Scatter plot of reconstructed versus generated values of W and \t\ for Monte Carlo 
events. 
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Figure 3: (a) The invariant 7r + 7r~ mass spectrum for all two track events selected by the 
untagged photoproduction trigger, (b) The distribution after applying all cuts, except those on 
W, Mn-,,- and pj>. The shaded area in both plots (hardly visible in (b)) is the mass distribution 
of the like-sign pion pairs, (c) The distribution after all cuts, including those on W and p\ 
(but not that on M^). 
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Figure 4: Acceptance as a function of (a), W (b) and (c) obtained with the DIPSI and 
the LEVEC generators. The quantity F plotted in (d) is defined in the text. Only statistical 
errors are shown. The horizontal bars indicate the size of the bins. 
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Figure 5: The energy spectrum in the forward calorimeter for: (a) elastic production (ep — ► 
ep°p) simulated with DIPSI; (b) inelastic production (ep —>■ ep°X) simulated with PYTHIA; (c) 
data; (d) a mixture of 89% elastic and 11% inelastic Monte Carlo events (histogram) compared 
with data (dots). 
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Figure 6: The acceptance corrected differential cross section da/dM^ for 60 < W < 80 GeV 
and \t\ < 0.5 GeV 2 . The points represent the ZEUS data and the curves indicate the results 
of the fit to the data using expressions fl7]-|T0"l). The dotted line shows the contribution of the 
Breit-Wigner term and the dash-dotted line that of the interference term; the shaded band 
indicates the size of the statistical uncertainty on the background term fps- The solid curve is 
the sum of these three terms. Only statistical errors are shown. The horizontal bars indicate 
the size of the bins. 
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Figure 7: (a) The integrated cross section cr 7P _^ p o p as a function of the centre-of-mass energy 
W. The ZEUS result (labelled "ZEUS 1993") is that given in flTJ) for the range 2M n < 
M W7T < M p + 5r , \t\ < 0.5 GeV 2 . The ZEUS result |2[ is also shown. The dashed line is the 



parametrisation from fl65 |. The vertical error bars of the ZEUS points indicate the quadratic 
sum of statistical and systematic errors. The horizontal bar indicates the size of the W region 
covered by the present measurement, (b) Same as (a) except that the results shown were 
obtained with the Spital and Yennie method. 
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Figure 8: (a) The differential cross section da/dt for ^p — > p°p in the mass range 2M,,- < < 
M p + 5r and 60 < W < 80 GeV. (b) The differential cross section da/dt for 7p — > p°p for 
60 < W < 80 GeV, as obtained by applying the Spital and Yennie method in each t bin. The 
continuous lines represent the results of the fit with the functional form (p0| ). Only statistical 
errors are shown. The horizontal bars indicate the size of the bins. 
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Figure 9: The slope b' t as a function of W as obtained from fits with function (p0|). The 
ZEUS result is that found by fitting to the data of Fig. [5]a. The vertical error bars indicate 
statistical errors only. The horizontal bar indicates the size of W region covered by the present 
measurement. 
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Figure 10: The differential cross sections da/dcos9h and da/dcph- The continuous lines rep- 
resent the results of the fits discussed in the text. Only statistical errors are shown. The 
horizontal bars indicate the size of the bins. 
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